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CALCULATIONS OF SOLAR PLASMA INTERACTION WITH MAGNETIC FIELD OF 


CURRENT LOOP AND TWO PARALLEL INFINITE LINES OF CURRENT* 
by Daniel J. McKinzie, Jr., and Anthony L. Julius, Jr. 

Lewis Research Center 

SUMMARY 

An investigation was made of the forces resulting from the interaction of the uncon- 
fined magnetic field emanating from a space vehicle with the undisturbed solar plasma 
(solar wind). The integrated theoretical drag per unit area and the acceleration of a cir- 
cular current loop exposed to the solar wind at zero degree angle of incidence were de- 
termined. Also, the integrated theoretical drag and lift per unit area, the restoring mo- 
ment per unit area per unit length of span were determined for the case of two parallel 
infinite lines of current at a 45° angle of incidence. The calculations of these quantities 
were based on the theoretical results reported by Sprieter and Briggs (ref. 4). 

For the case of the current loop, the drag per unit area varied from approximately 
-7 -4 

7. 2x10 to 3. 35x10 newton per square meter over a range of proton number densities 
6 8 

from 10° to 10° particles per cubic meter and a magnetic field parameter from 10” 

_ 2 

to 10 tesla. An examination of the particle fluxes from a major solar flare showed it 
to have a negligible effect on the drag when compared with that from the solar wind. Ac- 
celerations were calculated for a hypothetical superconducting coil (hollow torus) that 

had a permanent unconfined magnetic field. The maximum acceleration calculated was 

-4 

approximately 1. 37x10 meter per second squared. For the case of the two parallel in- 
finite lines of current at a 45° angle of incidence, the lift to drag ratio was approximately 
1 to 6. 8. 


INTRODUCTION 

The discussion of the impact of a steady stream of plasma emanating from the Sun 
(solar wind) and interacting with the geomagnetic field of the Earth led Chapman and 
Ferraro (ref. 1) to the concept of a hollow volume (magnetosphere) carved out in space 
from which the solar wind is excluded. The magnetosphere (named by T. Gold) may be 
defined as the region in the vicinity of the Earth in which the Earth’s magnetic field dom- 
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inates all particle dynamic processes. The boundary of the magnetosphere is called the 
magnetopause. In 1951 Martyn (ref. 2) described the steady plasma flow as exerting a 
constant pressure on the geomagnetic field whose own pressure maintains equilibrium. 

In 1958, Dungey (ref. 3) also pointed out that, for the case of steady-state conditions at 
the Earth’s magnetopause, the aerodynamic pressure of the deflected stream is balanced 
by the magnetic and electrical forces on the charged particles. 

Though the fundamental concept and equations governing the interaction between a 
neutral stream of ionized solar corpuscles and a three-dimensional magnetic dipole have 
been established for many years (ref. 4), the solution of the resulting equations has 
proved to be a difficult task. Most investigations of this problem have involved a reduc- 
tion in the number of dimensions, but Beard (ref. 5) has presented an analysis in which 
the full three-dimensional character of the problem is retained throughout. He was able 
to derive a nonlinear partial differential equation for the coordinates of the magnetpause. 
However, Beard's results for the meridional plane did not satisfy the governing differen- 
tial equation developed in the main body of the analysis. This fact was later pointed out 
by Sprieter and Briggs (ref. 4), and the problem was reworked for the meridional plane 
containing the dipole axis and the velocity vector of the incident stream for several cases 
in which these two directions were not perpendicular. The principal contribution of 
Sprieter and Briggs was the determination of the trace of the magnetospheric boundary in 
the meridional and equatorial planes for these several cases. The predicted profiles or 
limiting surfaces of the magnetosphere obtained in reference 4 are similar in scaled 
shape and size to the experimental results from reference 6 and also the satellite data of 
reference 7. The determination of the trace of the two dimensional magnetospheric 
boundary was also included in reference 4. 

While the analysis of Sprieter and Briggs for the three-dimensional case is based on 
the assumption that the Earth's magnetic field may be approximated by the magnetic field 
of a dipole, the present analysis is based on the assumption that the magnetic field of a 
dipole may be approximated by the magnetic field of a current loop at large distances 
compared with the radius of the current loop. This assumption applies to the two- 
dimensional analysis also. It can be shown that, for the three-dimensional case (ref. 4), 
the maximum downstream cross-sectional area of the magnetosphere may vary from 
several to hundreds of times the projected area enclosed by the current loop, depending 
on the magnitude of its magnetic moment and the momentum of the plasma stream. Thus, 
the momentum exchange at the magnetopause produced by a current loop that is part of a 
space vehicle could have a significant effect on the motion of that vehicle. The possible 
magnitude of this effect is investigated in this report. 

The theoretical analysis of reference 4 and its tabulated results are applied to de- 
termine the integrated momentum exchange over the magnetopause for two hypothetical 
space vehicle configurations. One configuration consists of a current loop (a three- 
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dimensional case), and the other consists of two parallel infinite lines of current (a two- 
dimensional case). The drag, lift, and moment are calculated for these cases. Finally, 
an estimate of the acceleration experienced by a nonlifting, hypothetical, hollow toroidal 
system is made. This system was studied for other purposes by Levy (ref. 8). All ap- 
plications made in this report are for the condition that the relative velocity between the 
vehicle and the solar wind is equal to that of the solar wind with respect to the Svui. 
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SYMBOLS 

area, m^ 

radius of current loop for three-dimensional case, m 

one-half normal separation distance for two-dimensional case, m 

total magnetic field vector, T 

permanent magnetic field vector, T 

three-dimensional magnetic field parameter, T 

two-dimensional magnetic field parameter, T 

total (tangential) magnetic field at magnetopause, T 
induced magnetic field vector, T 

unit length of span of two-dimensional magnetosphere (see fig. 3), m 

speed of light, 3x10^ m/sec 

total lift or drag force (see fig. 6), N 

total drag force (three dimensional), N 

total drag force (two dimensional) , N 

drag force (two dimensional) over upper magnetopause, N 
drag force (two dimensional) over lower magnetopause, N 
total lift force (two dimensional), N 

lift force (two dimensional) over upper magnetopause, N 
lift force (two dimensional) over lower magnetopause, N 
force (three dimensional) normal to magnetopause, N 
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force (two dimensional) normal to magnetopause, N 

2 

gravitational constant 9. 814 m/sec 
current, A 

arbitrary constant of integration 

total resultant moment acting about center of two-dimensional case 
(see fig. 4), N-m 

total moment resulting from drag forces about center of two- 
dimensional case (see fig. 4), N-m 

moment resulting from drag forces of two-dimensional case in cp = tt/2 
plane (see fig. 4), N-m 

moment resulting from drag forces of two-dimensional case in 
cp = -tt/2 plane (see fig. 4), N-m 

total moment resulting from lift forces acting about center of two- 
dimensional case (see fig. 4), N-m 

moment resulting from lift forces of two-dimensional case in cp - tt/2 
(see fig. 4), N-m 

moment resulting from lift forces of two-dimensional case in cp = -tt/2 
plane (see fig. 4), N-m 

2 2 

magnetic moment of circular -current loop, ttsl I, A-m 

-27 

mass of positive ion (proton), 1.67X10 kg 

number density of positive ions (protons) in undisturbed solar wind, 
particles/m 0 

unit vector in direction of outward normal to magnetopause 

- R 9 

solar light pressure at Earth’s orbit, 9x10 N/m 
coordinate of permanent magnetic field and magnetopause, m 
unit radial vector 

distance from apex of magnetosphere to center of three-dimensional 
current loop for special case in which A = 0, m 

distance from apex of magnetosphere to center of two parallel infinite 
lines of current for special case in which A = 0, m 

surface distance measured in meridional or equatorial plane, m 

solar constant (eq. (47)) 1. 35x10^ J/(sec)(m^) 
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stagnation point in two- and three-dimensional cases, m 

cutoff point of integration on three-dimensional magnetopause, where 
| cos ^ | sw 0. 03 (see fig. 1), m 

velocity of positive ions (protons) in undisturbed solar wind relative to body 
considered, m/sec 

velocity vector, m/sec 

coordinate of permanent magnetic field and magnetopause 
cutoff point of integration 

cutoff point of integration for lower magnetopause 
cutoff point of integration for upper magnetopause 


unit meridional vector 
9 - X sin <p (see fig. 4) 

value of 9 at pseudostagnation point (see figs. 3 and 4) 

angle between direction of incident stream and plane of current system (angle of 
incidence) 

\d - 77/2 | (see figs. 3 and 4) 

-7 2 

magnetic permeability in free space 477X10 (kg)(m)/C 

dimensionless radial coordinate (p = r/r ) of magnetopause 


dimensionless radial coordinate 


p 2 = r /( r o). 


of magnetopause 


coordinate of magnetopause (see fig. 2) 

angle between velocity vector of undisturbed particles of solar wind and outward 
normal to magnetopause at their point of impact (see fig. 1) 


THEORETICAL CONSIDERATIONS 

Presented herein are some of the more fundamental assumptions and equations on 
which the work of Sprieter and Briggs (ref. 4) is based. The assumptions are as fol- 
lows: (1) A steady-state, neutral ionized corpuscular plasma stream flows radially out- 
ward from the Sim. (2) The geomagnetic field is represented by the field of a magnetic 
dipole. This magnetic field has the effect of creating a magnetosphere bounded by a thin 
current sheath (the magnetopause) that has the property of confining the solar plasma to 
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the exterior and the geomagnetic field to the interior of the magnetosphere. (3) Interior 
to the magnetopause 


V X B = 0 

where B represents the magnetic field. (4) The particles of the solar stream move in 

straight lines to the magnetopause where they are, in effect, specularly reflected and 

and returned to the stream. In so doing, they exert an aerodynamic pressure of 
2 2 

2mnV cos i// on any element of the boundary for which cos ^ 0, where x p represents 
the angle between the velocity vector of the undisturbed particles of the solar wind and 
the outward unit normal vector n at the point of impact (see fig. 1) and m, n, and V 
represent the mass, number density, and velocity, respectively, of the positive ions rel- 
ative to the body considered in the undisturbed solar wind. (Fig. 1 is a modification of a 
figure in ref. 4.) (5) An element of the surface that fails to comply with the condition 
that cos s 0 is shielded from the stream and experiences no aerodynamic pressure. 

Dungey (ref. 3) showed that the aerodynamic pressure is balanced by the magnetic 
pressure / 2y u , where B is the total (tangential) magnetic field in teslas, at the 

magnetopause. This pressure balance determines the location in space of the magneto- 
pause. These considerations lead to the following relation that must be satisfied at the 
magnetopause. 


— — = 2mnV 2 cos^V (1) 

2 ^o 

(The equations developed in this report are written in the mks system of units. ) 

The total magnetic field B in the interior of the magnetosphere is the sum of the 
permanent magnetic field Bp and the induced magnetic field B' due to the currents in 
the magnetopause. As pointed out in reference 4, this field depends on the shape of the 
magnetopause, and its properties are described by the solution of the following magnetic- 
field equations 


div B = 0 

(2a) 

curl B = 0 

(2b) 


which satisfy the boundary conditions that the normal component of B vanish and the 
tangential component of B be given by equation (1) at the magnetopause. For the three- 
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2 2 

dimensional case (current loop), the permanent magnetic field where r » a is given 

by 



11 M 

P (§ sin 9 + ?2 cos 9) 

4ni 


(3) 


where represents the magnetic moment of a circular current loop, a represents the 
radius of the current loop, I represents the current in the loop, r represents the radial 
coordinate of the permanent magnetic field, and r and 9 represent the unit coordinate 
vectors. In this report, Bp Q is defined as 


M 

B S P 0 
P° 3, 
a 4ir 


(4) 


where B has units of magnetic field and is used as a parameter. As presented in ref- 
erence 4, B is a measurable quantity (i. e. , the magnetic field at the surface of the 
Earth), but herein, Bp Q as defined in equation (4), is simply a convenient parameter and 
does not necessarily specify the magnetic field at a or at any other specified point in 
space. Substituting into equation (3), yields 



(§ sin 9 + r 2 cos 9) 


(5) 


For the two-dimensional case, 


B 


P 


— — - (§ sin 9 + r cos 9) 
ra 2r 2 


( 6 ) 


where 2a„ represents the normal separation between the two parallel infinite lines of 
* 2 2 

current, and where the condition r » must be satisfied. The magnetic field pa- 
rameter for the two-dimensional case may be defined as 



where ii 0 l/4ir&2 is the magnetic field due to an infinitely long straight line of current 
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at a point of separation equal to 2ag. Substituting equation (7) into equation (6) results 
in 



( 8 ) 


It is assumed in reference 4 and in this report that the total tangential magnetic field at 

the magnetopause B is equal to twice the magnitude of B at the magnetopause. It 

s P 

should be noted that any magnetic system whose magnetic field may be expressed in 
either of the preceding forms (eqs. (5) or (8)) is similar to these systems. 

In reference 4, the following differential equations were obtained for the location of 
the magnetopause in the meridional plane containing the dipole axis and the Sun-Earth 
line: 


and 


dp _ I p ■ sin[0 - A sin(±7r/2)] - sin 0 
^ ' p^ cos[0 - A sin(±7r/2)] + 2 cos 9 


dp _ J p 6 sin [0 - A sin(±7r/2)] + sin 6 
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p'* cos[0 - A sin(±7r/2] - 2 cos 0. 


(9) 


(10) 


These equations differ only in the signs of the last term in the numerator and denomina- 
tor on their right sides. The solutions of these equations result in two families of 
curves that have been matched through physical reasoning to describe the complete pro- 
file of the magnetopause in the meridional plane. In equations (9) and (10), p is the di- 
mensionless radial coordinate (p = r/r Q ) of the magnetopause, A represents the angle 
between the direction of the incident stream and the normal to the axis of the Earth, and 
0 is a spatial coordinate of the permanent magnetic field. The angle ±n/2 represents 
the value of cp in the meridional plane (see fig. 2, a modification of a figure from 
ref. 4). The distance r Q is measured from the apex of the magnetosphere to the center 
of the dipole for the special case in which A = 0. For the three-dimensional case 
(ref. 4), 
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o p 


(47r)^mnV^ 


1/6 


b: 


1/6 


= a 


po 


{li Q mnV 


The solutions of equations (9) and (10) are given by 


p cosf# - A sin(±7r/2)l — — cos 9 = K 

p 2 


and 


( 11 ) 


( 12 ) 


p cos[0 - A sin(±7r/2)] + — cos 9 = K (13) 

P 2 

where K is an arbitrary constant of integration. Only one solution from each family 
matches properly in the meridional plane to give the physically continuous profile shown 
in figure 1. 

In reference 4, the final form of the differential equation for the shape of the mag- 
netopause in the equatorial plane (0 = n/2) is given. It is written for the special case 
where A = 0 and <p is restricted to the interval n/2 < <p < 3/2 n: 


dp 
d cp 



sin cp cos <£ + 


6 2 
p cos cp - 



(14) 


The solution of this equation gives the shape of the magnetopause in the equatorial plane 
(fig. 2). 

The foregoing is a presentation of some of the more fundamental assumptions and 

equations that made it possible for Sprieter and Briggs to describe the magnetopause of 

the Earth. They tabulated values representing the magnetopause as a function of 9 and 

2 2 2 2 

cp (see fig. 1). With the assumptions made herein that r »a and r » ag, the 
tabulated results of reference 4 are applicable to the cases of the current loop (three- 
dimensional case) and the two parallel infinite lines of current (two-dimensional case). 
The following is an analysis of the effect of the magnetic field solar wind interaction that 
gives rise to the forces over the magnetopause for the two- and three-dimensional cases. 
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Determination of Drag for Current Loop at Zero Angle of Incidence 


If specular reflection is assumed at the magnetopause, the force normal to an ele- 
ment of area dA of the magnetopause can be written 


dF n = 2mnV^(cos^)dA 


(15) 


Now, from figure 1, 


dA = 2i7r^(cos0)pd|— 


ft) 


(16) 


where s is the distance along the surface and is measured in either the meridional, or 
the equatorial plane. The drag force acting in the direction of the velocity vector V is 

dF D = -(cos i//)dF n (17) 


Substituting equations (15) and (16) into equation (17) and integrating yield 

cos \]s(cos9)pd 

where Sj is the cutoff point of integration on the magnetopause |cos xf/\ « 0. 03 (see 
fig. 1), (180 - i//) < 90°, 9 ^ 90°, and s q represents the pseudostagnation point s = 0. 
In polar coordinates, 




d 



+ 



d<9 


(19) 


When equation (19) is substituted into equation (18) and simplified, the results are 
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cos*V(cos 0) 


Hs)T 


<p=ir/2 


if 


cos ^(cos0)pp + 





( 20 ) 


where the upper limit 9 on the right side of equation (20) is the point at which the inte- 
grand becomes negligible and is the cutoff value 9 = (5/12) 7r given in reference 4 (see 

3 ^ 

fig. 1). The expression cos xf/ may be written as an explicit function of 9. The first 
integral on the right side of equation (20) is the contribution to the drag in the plane 
cp = jr/2 (fig. 1), while the second integral is the contribution in the plane cp = -n/2. 

In deriving equation (20), the magnetopause was assumed to be axisymmetric. Fig- 
ure 2 shows the shape of the magnetopause in the meridional and equatiorial planes. The 
figure indicates that the magnetopause is not axisymmetric, although the maximum down- 
stream cross-sectional radius perpendicular to the y/a-axis appears to approach an as- 
ymptote of equal magnitude in both the equatorial and the meridional planes. Because of 
this asymmetry, two separate calculations were made by using the meridional-plane 
profile cp = ±90° in one case, and the equatorial plane profile 9 = 90° in the other. 

The agreement between the two values of the integrals, obtained by using the trapezodial 
rule, was within ±1 percent. The mean value of the resulting integral was -0. 513. It is 
felt that there are no irregularities in the magnetopause in any other axial plane; thus, 
it is assumed that this value adequately represents the value of the integral. Substituting 
the mean value of the integral into equation (20), yields 

F d = 2. 05 mnV 2 7rr 2 (21) 


Determination of Lift and Drag for Two-Dimensional 
Case at 45° Angle of Incidence 

Sprieter and Briggs (ref. 4) determined the differential equations that describe the 
trace of the magnetopause in the meridional plane for the two-dimensional case. The 
equations are 
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(22) 


and 


dpg 

d& 



sin[0 - X sin(±7r/2)] - sin 9 
— 

sin [0 - X sin(±7r/2)] + cos 9 



p| sin [0 - X sin(±7r/2)] + sin 9 
0 

Pg cos[0 - X sin(±7r/2)] - cos 9 


(23) 


where pg is the dimensionless radial coordinate P 2 =r /( r 0 ) of the magnetopause, and 
( r o) represents the distance from the apex of the magnetopause to the centerline be- 
tween two parallel infinite lines of current for the special case in which X = 0. From 
reference 4, the expression for (r Q ) is given by 




(24) 


The solutions of equations (22) and (23) are given by 


Pg cos[0 - X sin(±7r/2)] - — cos 9 = K 

p 2 


(25) 


and 


P 2 cos[0 - X sin(±7r/2)] 


+ — cos 9 = K 

P 2 


(26) 


where K is an arbitrary constant of integration. Sprieter and Briggs (ref. 4) solved 
these equations for the magnetopause of the two-dimensional case and tabulated values 
representing the magnetopause as a function of 9 and P 2 in the meridional plane (see 
fig. 3). In the present analysis, only the tabulated values for X = 45° were used because 
of the larger anticipated lift forces for this case. As in the three-dimensional case, the 
change in momentum over the magnetopause of the two-dimensional case that incorpo- 
rates equation (19) is given by 
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I? 


d ( F n) 9 = 2mnb ( r o) 9 v2 COs2 ^ 


2 / dp 2\ 

P 2 + l 

\dO 


nV2 


de 


(27) 


where represents the local normal force acting on the differential area 




d<9 


of the surface of the magnetosphere, and b represents the length of the span of the two- 
dimensional magnetopause (see fig. 3). 

Lif t calculations. - The lift with respect to the wind axis (fig. 3) is given by 

d ( F L) 2 = sin ^ d ( F n) 2 < 28 ) 


where the positive direction of the lift force is indicated in figure 3. Substituting equa- 
tion (27) into equation (28) yields 


d 



= 2mnb 



V 2 cos 2 i^ sin xp 



(29) 


2 

where cos xf/ and sin xf/ are explicit functions of 0. Integrating equation (29) over the 
upper surface of the magnetopause shown in figure 3 results in 
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I 


where A g = 


- 7 r/ 2 1 and 9 is the value of 9 at the pseudostagnation point (fig. 3). 

_ .3 / S 

The upper limit n ) of the right side of equation (30) results in a negligible value of the 

' 'T ( \ 

integrand and is the cutoff value (3 j = 7r/6 given in reference 4. Integrating equa- 
tion (29) over the lower surface of the magnetopause yields 



where the lower limit (o ) of the right side of equation (31) results in a negligible value 

\ c/ B 

of the integrand and is the cutoff value (b \ = 77/3 given in reference 4. Combining 

\ c/ B 

equations (30) and (31) yields 


( f l ) 2 = ( f, 4 + ( f i \ 


(32) 


which results in the total lift force over the two-dimensional magnetopause. 

The evaluation of the integrals in equations (30) and (31) is 0. 2897. Substituting this 
value into equation (32) results in 



= 0. 5794 mnb 



V 


(33) 


Drag calculations , 
given by 


- The drag, in the direction of the velocity vector V (fig. 3). is 


d 




(34a) 
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where d(F n ) is given by equation (27). Substituting equation (27) into equation (34a) 
2 

and simplifying yield 


d 



-2mnb(r o ^ cos 
2 


p\ + 



(34b) 


Integrating equation (34b) over the upper surface of the magnetopause (see fig. 3) results 
in 



where 

( 9 c) t - 


For the lower surface 




(p= 7 t /2 


+ 



COS^lp 


p\ 


+ 





d6> 


j 


(p=-ir/2, 


(36) 
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where 


(e c ) B » »/s 

Combining equations (35) and (36) results in 

( F d ) 2 ' ( F b ) 2 + ( F b ) 2 < 37 > 

which represents the total drag force over the two-dimensional magnetopause. The 
evaluation of the integrals in equations (35) and (36) is -1. 976. Substituting this value 
into equation (37) results in 


(f d ) 2 = 3. 952 nrnb(r 0 )_V 2 


(38) 


Determination of Moment for Two-Dimensional 
Case at 45° Angle of Incidence 


The moment about the axis of the two-dimensional case is illustrated in figure 4, 
which is a modification of a figure from reference 4. The assumption made in calcu- 
lating this moment is that the two parallel infinite lines of current have been perturbed 
in such a manner that their angle of incidence is 45° to the velocity vector of the solar 
wind. In reference to figure 4, the moment resulting from the drag force is 



2(4*»/2) 


( m d). 


2(cp=-Tr/2) 


Now, in the plane <p = tt/2, 


(39) 



where 



In the plane cp = -tt/2, 


( m d) 


2{<p=-ir/2) 


- 2mn ( r ol bv2 


' /*( e c) T 


p 2 cos cos xf/ 


2 / dp 2\ 

Po + I ) 

V ddl 


n 1 / 2 


de 


f 

J(e„ 


+ i P2 cos e^ cos ^ 

H 

where = n/ 6 and (o ^ -tt/2. 

T B 

The integrated moment resulting from lift forces is 


2 l dp 2 


fll/2 


P 2 + 


de 


de 


(m t ) = (m t ) +(m t ) 

V L ) 2 \ ^/ 2 {cp=tt/2) V ^ J '2(<p=-ir/2) 


In the plane cp = tt/2, 


(m t \ = 2mn(r^) bV 2 

V h >2{cp=TT/2) V0/ 2 


r /• tt/2+\ c 

Jo 


P 2 sin 0 R sin if/ cos *xf/ 


2 M dp 2 

P2+ U. 


nV2 


de 


/%TT 

r / x ? 

1/2 

[ 2 

I p >2 sin sin \p cos \js 

2 dp 2\ 
p 2 + [ ) 
\ dd / 

de J 

Jtt/2+X 




In the plane <p- -tt/2, 



2 2 2 
p 0 sin sin if/ cos xf/ p„ + — - 

, V d0 , 


where (d ^ = tt/ 6 and (d = 77/3. 

Therefore, the total resultant moment is given by 


M 2 - ( M d ) 2 + ( M l ) 2 < 45 

Because 2mnbV^(r^ is a common factor in equations (40) to (45), all the integrals 

2 

may be summed. The value of their summation is approximately 0. 066. The resultant 
moment acts as a restoring moment (fig. 4). Thus, the total moment evaluated for the 
two-dimensional case at a 45° angle of incidence to the solar wind is given by 


M 2 » 0. 1320 mnb( 


DISCUSSION OF RESULTS 

The results of calculations made by using the preceding analyses are presented in 

figures 5 to 7. Figure 5 is a plot of the magnetic field parameter Bp Q against the ratio 

of the drag to a^ for the current loop at a zero degree angle of incidence to the solar 

wind. The curves are plotted for a constant proton number density n of the solar wind. 

6 7 8 

The number densities considered are 10 , 10 , and 10 particles per cubic meter. The 
number density is implicit in the relation between Bp Q and the drag per unit area. The 
data were computed by assuming the mean relative velocity of the undisturbed solar wind 

5 

to be 5x10 meters per second. An examination of the effect of a major solar flare 

showed that it makes little or no contribution to the magnetic interaction forces because 

of its small particle flux densities. Also shown in figure 5 is the value of the ratio of the 
2 

drag to a that results from light pressure acting on a circular normally oriented re- 
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flective surface having a radius a equal to that of the current loop. The resulting drag 
was calculated by assuming the value of the solar light pressure at the Earth's obrit to be 


P 


o 


— - = 9xl0" 6 N/m 2 
c 


(47) 


where 


S Q = 1.35X10 3 j/(sec)(m 2 ) 


and 


c = 3xl0 8 


m/sec 


The ratio of the drag per unit area (due to the solar wind magnetic field interaction) 
to that of the solar light pressure varies from a minimum of approximately 2. 5xi0 -2 to 
a maximum of approximately 12 over the range of the magnetic field and number density 
considered. Thus, if the magnetic moment of the current loop is large enough, the in- 
teraction of the solar wind with the magnetic field of the loop can result in a larger drag 
force than that due to solar light pressure. 

In order to assess the navigational possibilities of the forces resulting from the in- 
teraction, a two-dimensional case was analyzed because of its mathematical simplicity. 
Calculations were made for a 45° angle of incidence to the solar wind. These indicated 
that the two lines of current would experience a lift to drag ratio of approximately 
1 to 6. 8. The curves of lift and drag per unit area, respectively, are shown in figure 6 

and are plotted against (El ) . The unit area referred to in this case is the unit length of 

\ po/2 


span b multiplied by half the separation distance ag of the two infinite parallel lines of 
current (see fig. 3). As in the three-dimensional case, the curves are plotted for con- 
stant proton number density n. Because the lift and drag forces are a function of the 
size and shape of the magnetopause, a multiplicity of field sources could possibly be ar- 
ranged to create various desirable shapes of the resultant magnetopause. The overall 
shape may range from a thin high-lift wedge to a blunt high-drag hemisphere. 

The existence of lift and drag forces may imply a resultant moment. Figure 7 is a 

plot of (b ^ against the restoring moment per unit area per unit length of span b for 
\ P °/2 

the two-dimensional case. The unit area for this case is the square of one-half the sep- 
aration distance between the two lines of current. In calculating this moment, it is as- 
sumed that the two parallel infinite lines of current have been perturbed in such a manner 
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that their angle of incidence is 45° to the velocity vector of the solar wind. The curves 
are plotted for constant proton number density n. The existence of the restoring mo- 
ment suggests that a method of attitude stabilization of this configuration with respect to 
the solar wind is available since a moment only exists for an angle of incidence greater 
than zero. It also implies that an additional force is necessary to overcome this re- 
storing moment in order to utilize the lift in a navigational capacity. 

To gain an idea of the magnitude of the accelerations that might be experienced, a 

hypothetical toroidal system was chosen. Estimates of the mass and the radius a of the 

toroid were based on the hypothetical hollow toroid system of reference 8. The system 

presented in reference 8 was devised to provide a toroidally shaped radiation-free 

volume. This volume would be shielded from solar and galactic corpuscular radiation 

having energies equal to and less than a specified value. In the present application, the 

system was designed to shield a volume of 10^ cubic meters from proton energies of 

200 MeV or less. The calculations were made by assuming a zero angle of incidence. 

The resulting mass, current, radius of the current loop, and magnetic field parameter 

-4 

were approximately 3x10 kilograms, 4x10° amperes, 3. 3 kilometers, and 3. 80x10" 
tesla, respectively. 

The application of the analysis of the geomagnetic field interaction with the solar 
wind (ref. 4) to a vehicular system such as that of reference 8 is subject to the condition 
that similarity criteria be satisfied. These criteria are set forth in reference 9. Of con- 
cern in the simulation is the condition that no ring currents exist (V x B = 0) within the 
magnetosphere of the hypothetical toroidal system, and also that the ion cyclotron radius 
of the protons at the magnetopause be less than the standoff distance of the magnetopause. 
The second, ion cyclotron radius, requirement is very important because it is at the 
magnetopause that the exchange of momentum from the solar wind takes place. These 
requirements have been satisfied in the calculation of the accelerations shown in table I: 


TABLE I. - THEORETICAL ACCELERATIONS OF TOROIDAL SYSTEM 


Three-dimensional 
magnetic field 
parameter, 

B po’ 

T 

Density, 

n, 

protons/m 3 

Acceleration of 
toroidal system a 

Sunlight reflective 
acceleration system b 

m/sec 3 

g 

m/sec 3 

g 

3.80X10" 4 

5xl0 6 

1. 85xl0' 5 

1.88X10' 6 

3X10" 3 

3.06X10' 4 

3.80 

10 7 

2. 95 

3.01 

1 3 

3.06 

3. 80 

t-* 

o 

OO 

1. 37xl0' 4 

1.40x10' 5 

3 

3.06 


a Ref. 8. 
^Ref. 10. 
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The accelerations are compared with the acceleration presented in reference 10 which 
results from the reflection of sunlight at 1 astronomical unit by a sheet of material with 
a specific gravity of 1.18 and a thickness of approximately 7.61x10” meter. The area of 
the sheet is equal to that of a disk which has the same maximum radius as that of the to- 
rodial system, and no mass has been included for the structural support of the sheet. 

The maximum calculated value of the acceleration of the hypothetical toroid considered 

-4 -5 1 

is 1. 37x10 meter per second squared (~10 g). This value is lg orders of magnitude 

lower than that experienced by the sheet in sunlight at 1 astronomical unit. 


SUMMARY OF RESULTS 

The results obtained from an investigation of the forces resulting from the inter- 
action of a magnetic field possessed by a hypothetical space vehicle with solar wind are 
as follows: 

1. The ratio of the drag due to the interaction of the magnetic field of a current loop 

with the undisturbed solar wind to solar light pressure drag imparted to a reflective sur- 

_2 

face that has an area equal to the area of the loop is approximately 2. 5x10 to 12 over 
the range of proton number densities and magnetic field considered. 

2. For the two-dimensional case at a 45° angle of incidence, the ratio of the lift to 
the drag is approximately 1 to 6. 8 and a restoring moment exists. 

3. The maximum calculated value of acceleration of the hypothetical toroid con- 

-4 -5 1 

sidered is 1. 37x10 meter per second squared (~10 g). This value is l-^ orders 

of magnitude lower than that experienced in sunlight at 1 astronomical unit by a sheet of 
very light reflective material for which no structural supporting weight has been in- 
cluded. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, February 3, 1967, 

129-01-05-05-22. 
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Figure 1. - Three-dimensional magnetosphere showing cross- 
sectional view in meridian plane at zero degree angle of 
incidence. 



igure 2. - View of coordinate system and boundary of equatorial plane and meridian plane 

containing current loop axis and Sun-Earth line at zero degree angle of 

incidence. 






Magnetic field parameter, 



Figure 4. - Cross-sectional view showing meridian plane and moment center for 
two-dimensional case at 45° angle of incidence. 



Figure 5. - Magnetic field parameter plotted against ratio of drag 
to current-loop radius at zero degree angle of incidence for 
current loop. 




Figure 6. - Magnetic field parameter plotted against ratios of 
lift and drag per unit area per unit length of span at 45° 
angle of incidence for two parallel infinite lines of cur- 
rent. 



Figure 7. - Magnetic field parameter plotted against restoring 
moment per unit area per unit length of span at 45° angle 
of incidence for two parallel infinite lines of current. 
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